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Abstract 

In this paper we study the behavior of (laser-cooled) m-atoms trapped in a 
cavity interacting with a photon • • • Cavity QED • • • and attempt to solve the 
Schrodinger equation of this model in the strong coupling regime. In the case of m 
= 2 we construct Bell-Schrodinger cat states (in our terminology) and obtain with 
such bases some unitary transformations by making use of the rotating wave approx- 
imation under new resonance conditions containing the Bessel functions, which will 
be applied to construct important quantum logic gates in Quantum Computation. 
Moreover we propose in the case of m = 3 a crucial problem to solve on Quantum 
Computation. 

This paper is a continuation of ^^1; CHI- We consider a unified model of the interaction 

of the two-level atoms and both the single radiation mode and external field (periodic 

usually) in a cavity. We deal with the external field as a classical one. As a general 
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introduction to this topic in Quantum Optics see jT], j2], 0. Our model is deeply related 
to the one of trapped ions in a cavity with the photon interaction (Cavity QED). 

In our model we are especially interested in the strong coupling regime, [Hj, JB], [T^ - 
One of motivations is a recent very interesting experiment, pi]. See |3] and [H] as a 
general introduction. 

In jH] and ^Hj we have treated the strong coupling regime of the interaction model of 
the two-level atom and the single radiation mode, and have given some explicit solutions 
under the resonance conditions and rotating wave approximations. 

On the other hand we want to add some external field (like Laser one) to the above 
model which will make the model more realistic (for example in Quantum Computation). 
Therefore we propose the unified model. 

We would like to solve our model in the strong coupling regime. Especially we want to 
show the existence of Rabi oscillations in this regime because the real purpose of a series 
of study ([H], [E], [IS], [in]) is an apphcation to Quantum Computation (see [T2j as a 
brief introduction to it). 

Since in the previous paper ^S] we have treated one-atom case, so we deal with two- 
atoms case in this paper. To solve the Schrodinger equation in this regime we construct (in 
our terminology) Bell-Schrodinger cat states, and make use of rotating wave approxima- 
tion under new resonance conditions containing the Bessel functions and etc, and obtain 
unitary transformations which are necessary to perform quantum logic gates. 

Our solutions might give a new insight into Quantum Optics or Condensed Matter 
Physics as well as Quantum Computation. 



We here make some preparations for the latter. Let {o"i, o"2j cts} be Pauli matrices and 
I2 a unit matrix : 

/ 



0-1 



1 

1 



0-2 




0"3 




1 

1 
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and cr+ = (l/2)((Ti + i<72), <J- = (l/2)(cri — 10-2). Let W be the Walsh-Hadamard matrix 



W 



V2 



1 1 
1 -1 



(2) 



then we can diagonahze ai as di = Wa^W ^ = Wa^W by making use of this W. The 
eigenvalues of di is { 1 , — 1 } with eigenvectors 



|1) 



V2 



-1) 



J_ ( 1 

V2 



(3) 



We have treated one atom with two-level, so we would like to generalize the method 
developped in to m atoms with two-level interacting both the single radiation 

mode and (classical) external fields like (m atoms trapped in a cavity) 



|0) |0) 



Then the unified Hamiltonian can be given as 



A 



i/i = wlAf®L3+^l5]dr^®(L++L_) + -5](T^^^®lL+^2EcOs(c^jt + </)j)(^l (4) 

j=l ^ j=l j=l 



where M = and a^^^ {k = 1, 3) is 



n^^'> - 1 
CTj. — U 



I2 <S> cTfe ® h 



I2 (j — position), 



see Here we are treating the following three cases at the same time : The notation 
{L+,L_,L3} means respectively 

' (N) {at,a,iV}, 
{L+,L„,L3}= <! (K) {K+,K_,Ks}, (5) 

(J) {J+,J-,J3}, 



and {a^,a,N} is the generators of Heisenberg algebra, {K+, K_, K^} and {J+, J_, J3} 
are a set of generators of unitary representations of Lie algebras 5^(1, 1) and su(2). 
They are usually constructed by making use of two harmonic oscillators (two-photons) 
{ai,a\}, {a2,4} as 



su 



1, 1) : = ai''a2\ K_ = 0201, = ^ (^ai^ai + a2"''a2 + l) , 

su{2) : J+ = ai^a2, J- = 02^01, "^3 ~ 2 ('^i^'^i ~ ^2^02) • 

See for example [13], [Ej in detail. 

We are interested in the strong coupling regime {gi ^ A), so we solve the Hamiltonian 
as follows (see [IS]). Let us transform (jU) into 



A 



Hl = ujIm ^L3 + J2 ^1 ® {9i{L+ + L_) + g2Cos{ujjt + 0^)1^} + ^ E ^3 ® 1l 



A 



a 



iJ) 



(6) 



where we have written Ho instead of Hlq for simplicity. First we diagonalize Hq. For that 
we set 

W = W ®W ®---^W G U{M) 
for in p. Then it is not difficuh to see (see [IB], [E]) 



= (W ® 1l) 



1m ® UJL3 + a'i^ (g) {gi{L+ + L_) + g2Cos{ujjt + 0j)1l} 



1m ® CUL3 + giJ2 ^3 ® (i^+ + i^-) + ^2 E ^3 ® cos{ujjt + (j)j)lL 



(w-i ® U) 

® 11) 



I r/£ I 

) ^ fiLs + (72 E A,cos(a;,t + 0,)U [ ((A| ® e^(^+-^-)) 



(7) 



where we have used the following 
Key Formulas 



(AT) uja^a + g^A{a^ + a) = Q e'^^'^'"'^) (n - j e^^-^^"") 

where fl = u, x = 2gi/u, 
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where Q = uo^l - (2^iA/cu)2, x = (l/A)tanh-^(2^iA/u;), (9) 
(J) cuJa + (7iA(J+ + J_) = 1] e-'^(-'+-^-) Js e'^(^+--'-) 



where = u^l + (2^iA/cj)2, x = (l/A)tair^(2^iA/c<;), (10) 
and have used concise notations A = (Ai, ■ ■ ■ , Aj_i, Aj, Aj+i, ■ ■ ■ , Am) and 



E= E E ■■■ E ; |A) = |Ai)®|A2)®---®|Am) ; A = A(A) = 5:A 

A Ai=±lA2=±l A„=±l i=l 



m 

\ 



for Aj = ±1. We leave the proof of Key Formulas to the readers. 

That is, we could diagonalize the Hamiltonian Hq. Its eigenvalues {En(t)} and eigen- 
vectors {|{A,n})} are given respectively 

(K(t), \{X,n})) = 

(N) ^](-^A2 + n) + (72E™lA,cos(cu,^ + 0,)> |A) ® e'^^^'-'^V), 
(K) r](ir + n) + (72Ef=iA,cos(cu,t + 0,), |A)®e-^(^+-^-)|i^,n), (11) 
_ (J) Q{-J + n) + g2ET=i^jCos{ujt + (Pj), |A) ® e-^(^+--'-)| J, n) 

for n G N U {0}, where En{t) = En{X) + g2 J2jLi XjCos{ujjt + (pj). Then Hq above can be 
written as 

^o = EE^"WI{A,ri})({A,n}| 

A n 

= EE I K(A) + (?2E Vos(^.-^ + 0,) \ |{A,n})({A,n}|. 

A n [ j=l J 

Up to this stage we could treat L = (N), (K), (J) at the same time. However in the 
following it is very difficult to treat them at the same time, so we shall focus our attention 
to the simplest case L = (N) ^. We leave the remaining cases as a future (challenging) 
task. 

Next we would like to solve the following Schrodinger equation : 
d / A 1 



^The main difficulty for L = (K) and (J) is the A dependence of x, see © and (fTn| . 



To solve this equation we appeal to the method of constant variation. First let us solve 
i-^^ — Hq^, which general solution is given by ^{t) — Uo{t)^o, where is a constant 
state and 



\{\,n}){{\,n}\. (13) 



A n 



The method of constant variation goes as follows. Changing like — ^o(^), we have 

J A / m \ A 

ij^^o = -^Uo^ {J2(T^^ ® ^Lj Uo^o = ^Hf^o. (14) 
Then making use of the fact (Ajjcra = (— Aj| and after long calculations, Hp becomes 

m 

jj^ _ ^ ^ y^gi|ta;(n-n')+ta;j:^(l-AjA(A))+2Aj(g2M)sm(a)jt+(/>j)| ^ 
A n,n' j=l 

l^n\e^Ma'-)\n')\{\n}){{\n'}^^^\, (15) 
where we have used the compact notations A(j) = (Ai, • • • , Aj_i, — Aj, A^+i, • • • , \jn) and 

(A(j)l = (All ® ■ ■■® (-Ajl ® ■■ ■® {\m\i 

= (A(,)| (8) (n|exp |^MM^(at _ a)| ^ 

and 

A(A(j)) = Ai H h Aj_i - Aj- + Aj+i H h A^ = A(A) - 2\j 

and 

K(A) - E^,{\^)) = |-^A(A)2 + n| - |-^A(A(,))^ + n j 
= u;(n - n') + a;x^(l - AjA(A)). 



In the following we set for simplicity — and 



em^g'^^^ for l<j<m, (16) 



then 

m 

^ ^ ^ ^eHM---')W(l-A,A(A))+2A,0,W}^^|gA,.(at-a)|^/^||^^^|^^|^^^/|^.^|_ (^7) 
A n,n' j=l 



Here we divide Hp into two parts Hp = Hp + Hp where 

m 

Hp = 5]5]X]e'{*""'(^-^^^(^»+'^^®^W}(n|e^^-("'-'^V)l{A,^})({A,n (18) 
A " i=i 

and 

m 

= E E Ee^^*"("-"')+*"^-'(^-^^-^(^))+'^^®^(*)}(n|e^^^^('^'"^^ 

A n,n' i=l 

(19) 

Noting (n|e'*'^"^""^|n) = (n|e~'^'^"^~"^|n) (Aj = ±1) by the results in section 3 of [T3] . 
Hp can be written as 

Hp =Y,{n\e^^'''-''^\n)\j2fl(i"^'''"'^'~^'''^''^^''^^^^^^^ \ ■ (20) 

n [ \ j=l J 

For us it is not easy to expand the equation above, so we consider the special case. 

Case of m=2 : 

2 

EE 

A j=l 

e^^"*"^'+'®^(*)^|{l,l,n})({l,l,n}(i)| + e^^-*--'+2e.W}||^ 



e^^*^^'+^«^^*^^|{l,-l,n})({l,-l,n}(i)| + e^i*-^'^-^«^^*^i|{l,-l,n})( 

e^^*^^"-'®^(*)>|{-ia,ri})({-l,l,n}(i)| + e^^*'^^'+2«^W>|{-l^ 

e^^-*"^'-'®^(*)^|{-l,-l,n})({-l,-l,n}(i)|+e*^-*--'-2®^^^^^ 
e^^~*-^-'+2®^W>|{l,l,r2})({l,l,n}(i)|+e*^-*-^''+2e2W}||^ 

e^^-*"^'-'®^(*)^|{-l,-l,n})({-l,-l,n}(2)|+e'^-*-^'-^®^W^ 

e'^*"^'-'^^(*)^|{l,-l,T2})({l,-l,n}(2)|+e^^*-^'+2«^W>|{l,-l,n^^ 

e^^*^^"-'®^(*)>|{-l,l,n})({-l,l,n}(i)|+e^^*-^'+2«^W>|{-l,l,n})({-l,l^ 
e'^"*'""'+2ei(t)}||i^l^^|^^_l^l^^|| ^gi{-t^x2+2e2(i)}||i^l^^|^^|l^_l^^|| ^ 

e*^-*""'~202{t)}||_i^ -1, n}){{-l, 1, n}| + e*{-*--'-2ei{t)}||_i^ n}){{l, -1, n}| + 

e*^*"^''-2®i W> |{-1, 1, n}) ({1, 1, n}| + e*^*-^''-2Q^W> |{1, -1, n}) ({1, 1, n}| + 
g.{W+2e.W}||_^^^^^|^^_^^_^^^||^g,{i..x2+2e,(t)}||^^_^^^|^^_^^_^^ (21) 



where let us once more remind 

|{l,l,n}) = |l)®e-"('^'-'^V), 
|{l,-l,n}) = |l)®|-l)®|n), 
|{-l,l,n}) = |-l)®|l)®|n), 
|{-l,-l,n}) = 1-1) ® 1-1) ®e^(''^-''V)- (22) 

Here we want to solve the equation i{d/dt)'^Q = Hp "^o completely, however it is very 
hard (see for example [Zj, [HI!, |2I|, |22j)- Therefore let us appeal to a perturbation theory. 
Now we use the well-known formula 

where Ja(x) are the Bessel functions. For a further simplicity we set 2g2/ujj = Tj. 

Here we define a kind of extended Bell states from ()22j) (|1) < — > |0), |— 1) < — > |1) in 
a usual manner) 

|{<^>i,^}) = ^(|{l,l,^}) + |{-l,-l,n})), 
I{<f2,^}) = ^(|{l,l,r^})-|{-l,-l,n})), 
|{<^>3,^}) = ^(|{-l,l,r^}) + |{l,-l,n})), 

|{*4,n}) = -^(|{-l,l,n})-|{l,-l,n})). (24) 

A comment is in order. We want to call these states (bases) Bell-Schrodinger cat states. 
Conversely we have 

l{l,l,^}) = ^(|{<^>l,^}) + l{<^>2,^})), 

\{-l,-l,n}) = ^m,,n})-\{<!>,,n})), 
\{-lAM) = ^ms,n}) + \{^,,n})), 

\{h-lM) = ^m3,n})-\{<^,,n})). (25) 



From (HH) 



However we are not interested in a transition from one state to an another one from 
the lesson obtaining one qubit case in ^^1; so we can remove Hp (n 7^ n) from the 
Hamiltonian. Next let us transform Hp. From ()21|) and ()23|) if we set 

A i=i J 
then we obtain after some algebras 

i^oF =(Jo(ri) + Jo(r2)) {e-^*"^'|{$i,n})({<|.3,n}| +e^*-^"|{$3,^ 

( Jo(ri) - Jo(r2) ) {e-^*-^' I {$2, ri}) ({<f 4, n} \ + e^'^^' \ {$4, n}) ({$2, n} \ } , (26) 

and 



e 



07^0 I 2 2 J 



e 



e 



a^Q I 2 2 J 

g^t^x^ ^ j Ja(-ri)- J^(ri) ^,,,,^ ^ J.(-r2)- J^(r2) ^,,,^ 1 |{$3,n})({<i>2,n}| + 

a^O I 2 2 J 

where we have used the basic relations 



Therefore 
d 



= Y^F^o = ^[-{n\e<'''--^\n)KoF + -(n|e-('^'-»V)^iF } ^o- (28) 

In the following we treat ^{n\e^^°'''~°'^\n)KQF a unperturbed Hamiltonian and the remain- 
ing a perturbed one. If we set for simplicity 

A 



2 
A 



He-(«'-")|n) (Jo(ri) + Jo(r2)) , 



EA,n,- = ^(r^|e"(«'-"V) (^o(ri) - Jo(r2)) , 



(29) 



then it is easy to solve 







= E[^A,n,+ {e-*--'|{$i,n})({$3,n}|+e^*'^-'|{$3,n})({$i,n}|} + 

n 

^A,n,-{e-^*'^^'|{$2,n})({$4,n}| +e^*--'|{$4,n})({$2,n}|}] ^o, 

see Appendix. 

As an application to Quantum Computation it is sufficient for us to consider one excited 
state. By making use of the method of constant variation again we can set ^'o(^) as 

*o(^) ^{Un,llC„,l{t) + Un,13Cn,3{'t))\{^l,'>T'}) + {Un,22Cn,2{t) + Un,24Cn,4{'t))\{^2, n}) + 

(Mn,3lCn,l(t) + Mn,33Cn,3(^)) | {$3, «}) + (Mn,42Cn,2(t) + «n,44Cn,4(t)) | {$4, ?^}) (30) 



for some fixed n, where 

Unit) 



Vn{t) 



\ «n,31 Mn,33 / 
^ Un,22 ^71,24 ^ 



\ Un,42 Un,44: J 

from the appendix. 



Atwx 



Atwx 





f 

—it 


f ° 


EA,n,+ 




> 


exp < 






> 






\ EA,n,+ 




) 






f 

—it 


( " 


EA,n- 




> 


exp < 






> 








UJX^ 


I 


> 



(31) 



(32) 
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Then substituting (plj) into and some tedious calculations lead us to 
( . \ ( 



with 
A. 

B 

C ■ 



.d 

I— 
dt 



c„,i 

Cn,3 
Cn,2 



Wn,31 "^71,33 



^ A ^ 



A C 
C D 
S 



1tn,ll Mn,13 
Wn,31 Wn,33 



Cn,3 
Cn,2 



(33) 



Ja(ri) + j„(-ri 



^gitaa;i _|_ 



e-'*'""'v4o, 
9 A 

2 ^ 

e-'*""'So, 



-itojx^;^ / I x(at-a) 
2 ^ ' 



= e 



D = e 



-itujx 



-itiox 



Co, 

2A , 



(n e 



x{at — a) 



itauji 







-r2) 




2 






-Ja{ 


-r2) 




2 




Ja{T2) 


-Ja{ 


-r2) 




2 




Ja{T2) 




-r2) 



itOLL02 



itauj2 



itauj2 



itauj2 



e-^*^^'Do. 



(34) 



We note that Aq = Aq, Dq = Dq and Bq = —Bq, Cq = —Cq (due to the fact J-a{x) 
Ja{~x)). 

The above matrix equation decomposes into two parts 



d 
I— 
dt 



Cn,l 
V Cn,3 



Un,31 Un,33 



-1 



\ Un,31 Un,33 



A 

A 

B 

C 



Un,31 Un,33 



Un,22 W„_24 
\ Mn,42 ^in,44 



Cn,l 
Cn,3 

Cn,2 
Cn,4 



+ 
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By the way, from (j!?T|) . (jH^ and 



exp < 




exp < it 



Doexp < it 




^A, 



E 



A,n, 



LUX 



Un,ll Un,13 
Un,31 Un,33 

UnA2 UnM 



I 



CnA 



+ 



\ ^n,3 

( r 

(^n,2 




exp < —it 



exp < 



-it 












exp < 














;) 


exp < 


!-( 










LOX 




+ 



(35) 



+ 



(36) 



At this stage we can set several resonance conditions and obtain corresponding equa- 
tions and solutions by making use of the rotating wave approximation like in ^21, |IS]- 
For example, let us consider the term 
/ 

exp I it 





y -EA,n,+ 

then from ()57j) in Appendix 




5o . 

exp i-it 
Co j I V EA,n, 












( 


Bo ' 
















1 Pn,- 








v 






^ Co 







v 
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Here for example, 



gi*(Mn,+ -/in,-)p 



ijt(t'n,+ -/in, + ) 



p-1 







Bo ] 










1 


[ 


Co 







-it{Vn,--Hn,-) 



for 



Bn = — (nle 



^(at — a) 



2 
A 



2 j' 



|n)E 



Now we set the resonance condition 

aU2 + Hn,+ - fJ'n- = <^ {-a)u;2 + Hn- - fJ'n,+ = 



(37) 



for some a ^ 0, then we can neglect the remaining terms by the rotating wave approxi- 
mation. That is, we obtain the time independent matrix 



A 







u r 












(' ) 


V 0. 




V 1 0. 







so we have 



d 



'dt 



(n e 



n) 



X 





(l 




p-1 








(l ) 




^ Cn,2 


\ 












p,..- 








\ 


oy 






0) 




[ ") 




\ C„,4 


/ 



^di 



\ Cn,4 



A 

2" 



— (n|e '\n) — — 



X 



p-1 







^ c ,\ 



where we have used the relation 

J_„(r2) - J-a(-r2) MT2) - U-r2) 
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After some algebras we obtain 




2 ^ ' ' ^ 2 



X 




,n,— fJ'n,+ 
(^A,n,+ + /^n,+)(^i,n,-+/^n, 

A, , „^, , Ja(r2) - Ja(-r2 



EA,n,+EA,n,- -E'A,n,+/^n,- 



-y(n|e^('^'-"V)- 



+ Ea ,n,— 
iEA,n,+ + l^l,+)iEA,n- + /^n ,-) 

from the appendix. 

If we prepare the vector notations 



X 



EA,n,+EA,n,- l^n,+EA,n,- 
\ EA,n,+l^n,- l^n,+l^n- 





' Cn,l 1 




' Cn,2 1 


Cod = 

















then the equation above can be written as 

.d 



'dt 



' Cod \ 




f 1 


f 


Cod ^ 


\ Cev J 








Cev J 



with 



\ Cn,4 



^ /^Ux(at-a)i \ ■^(^2) - Ja{-^2) EA,n,+f^n,^ + EA,n-f^n,+ 



and 



,n,— A'""., 

The (formal) solution is easily obtained to become 

/ 



EA,n,+EA,n,- EA,n,+f^n,- 
Ea ,n,— l-^n,+l-^n,— 



Cod{t) \ I int 

exp 

Cevi^^ 
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Next let us calculate the unitary transformation of ()4H). It is easy to see 

/ 



K 



1 



(^i,n,+ + /^n,+)(^A,n,-+/^n,- 



SO 



F2 



1 I -E'A,n.- I / \ 
2 (,^A,n,-, /Xn,-J 

1 (-^A,n - + A^n -) (-^A.n,-, /^n,-), 



P2 I ,,2 



-^A,n,-t 



Namely, K^K and KK^ are just projection matrices and moreover satisfy 

[KK^Y = KK^, [K'^KY = K^K, KK^K = K, K^KK^ = 

for n>l. Therefore after long algebras using the relations above we have 

K 







/ 




iTZt 


exp < 








[ ' 




(l 


- KK'^ + cos 


V 


isin {Jy) 



zsm 



(f)^ ^ 
1 - K^K + cos (f ) K^K J 



so the solution is explicitly 



/ 



Cod{t) 



1 -KK^ + COS (^]KK^ 

T 



\ Cev{t) 



in (f ) K 



isin (Jf) K 



(42) 



(43) 



(44) 



1 - K^K + cos (f ) K'^K 



Cod(O) 
Cet,(0) 



(45) 
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with (jig) and 

This is one of unitary transformations that we are looking for in Quantum Computation. 
We can also obtain another ones by setting different resonance conditions like (pTTj) (we 
leave them to the readers). 



Lastly let us summarize our result. By (0} the Hamiltonian becomes 

= E 1^^?^ + 92Cos{ujjt + 0j)ap^| ® 1l 



(46) 



if there is no photon interaction, so the dynamics of each qubit space is independently 
determined by 

Hlj = + 92Cos{ujjt + 0,)ap'^ . (47) 

Therefore the total space of m-qubits is just 



C^® •■■OC^® •■■OC^ where = Vectc{|0), |1)}, 

and unitary transformations of each qubit space are obtained by manipulating laser fields. 
However to solve the Schrodinger equation (neglecting the suffix) 

d . (A 



^di^ = |y^3 + 92Cos{ujt + 0)cri| ^ 



is not so easy, see for example [TUj . 

When m = 2, the interaction (driving by the photon) between two qubits is given 
by unitary transformations like ()45|) ■ ■ • controlled unitary gates including the controlled 
NOT. 
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This is just our scenario of Quantum Computation in the strong couphng regime. 



Case of ni=3 : 

We present a very important 

Problem : Let us consider three atoms in a cavity. How can we construct C-NOT (or 
C-Unitary) operations for any two atoms among them ? 
See the following pictures : 



C-NOT 




These constructions are very crucial in realizing quantum logic gates, however we 
have not seen such constructions in any references. We will attack this problem in a 
forthcoming paper. 

In this paper we treated the two-atoms case in a cavity QED and constructed unitary 
transformations by making use of the rotating wave approximation under new resonance 
conditions containing the Bessel functions. 

17 



These will be applied to construct several quantum logic gates in Quantum Compu- 
tation. Moreover we would like to treat a general case, which is at the present not easy 
due to some technical reasons. 

By the way, according to increase of the number of atoms (we are expecting at least 
m = 100 in the realistic quantum computation) we meet a very severe problem called 
Decoherence, see for example [Oj and its references. We don't know how to control this. 
One way protecting against this may be to deal with N-level system (then we can reduce 
the number of atoms in a cavity). A generalization of the model to N-level system (see 
for example jHj, jTHI, [201; |21I) is now under consideration and will be published in a 
separate paper. 
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Appendix : Some Useful Formulas 



In this appendix we solve the following equation 



= aHijj, 

at 



(48) 



where a is a constant and 



= e-*'^*|l)(-l|+e*'*|-l)(l| and ^ = a{t)\l) + b{t)\-l) . 



(49) 



Then it is easy to get a matrix equation on {a, b} 

I 



d I a 



V 



ae-^^* ] 






ae^^* J 




: 



d ~ 



(50) 
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Noting the decomposition 



ae 



ae 
m Q 



iet 



V 


(l 


\ 


(o 






\ 








^ a 


J 


\ 





then the equation becomes 



.d J 



a 

a e 



\ 



ijj for ijj 



-iet 



The solution is easily obtained to become 

/ 

V 

where (oq, 6o)^ is a constant vector and 




fl 


\ 


exp < 


-it 


U 




1 


\ 








^ a 







d 

i—U = HU. 
dt 



If we set 



Uu Ui2 



\ 



U21 U22 J 

(2 corresponds to —1) then above can be written as 

= {uiiao + Ui2bo) 1 1) + {u2iao + ii22&o) I -1) 

with constants {qq, bo}- 

In the method of constant variation in the text we change like Qq — > cio{t) 
bo{t). 

Let us make some comments. For 



a 6 J 



we can easily diagonalize A as follows : 

^ a a ^ 



V 



y/a^+lJ? Va2+i.2 
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where 



Therefore we obtain 



Q{t) = e-'*^ = 



JL 



For the simphcity we set 



-itfj, 



\ ^a2+^2 Va2+i.2 y 



e 0(2). 



In the text this is used as 



/ -BA,n,( 



for cr = ± 



with 



for 7 = Then 












( 


On 


ai2 


\ 






;) 












V 


1 






\ 


021 


022 


/ 



where 



On 



021 



,n,+IJ'n,- ,n,—l^n,-\ 







+ \jEA,n, 


- + /^n,- 






+ EA,n,- 





ai2 



022 



EA,n,+^n,- + Ea 



EA,n,+ l^n, 



Eln,- + <■ 



EA,n,+^n- + EA,n-^n, 



(1„ 
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